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Abstract 



""O The calorimeter, range detector and active target elements of the T2K near detectors rely on the Hamamatsu Photonics Multi- 

c/5 Pixel Photon Counters (MPPCs) to detect scintillation light produced by charged particles. Detailed measurements of the MPPC 

Ch gain, afterpulsing, crosstalk, dark noise, and photon detection efficiency for low light levels are reported. In order to account for the 

^ impact of the MPPC behavior on T2K physics observables, a simulation program has been developed based on these measurements. 

Q The simulation is used to predict the energy resolution of the detector. 
• i— I ■ 

^ 'Keywords: photosensors, photodetectors, multi-pixel avalanche photodiodes, scintillator 
J^ PACS: 29.40.Mc, 29.40.Wk, 29.40.Wj 
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1. Introduction 

The Tokai to Kamioka (T2K) project [1] is a second- 
generation long-baseline neutrino oscillation experiment that 
uses a high intensity off-axis neutrino beam produced by the 
30 GeV proton beam at the Japan Proton Accelerator Research 
Complex (J-PARC). The first phase of the T2K experiment pur- 
sues two main goals: a sensitive measurement of On, and deter- 
mination of the parameters sin 2 2$23 and Am 2 3 to better accuracy 
than any previous experiment. 

To reach these physics goals, precise knowledge of the neu- 
trino beam flux and spectrum, and the neutrino interaction cross 
sections is required. To perform the required measurements, the 
near detector complex (ND280 [2]) was built at a distance of 
280 m from the hadron production target. The complex has two 
detectors (Fig. 1): an on-axis detector (neutrino beam monitor), 
and an off-axis neutrino detector located along the line between 
the average pion decay point and the Super-Kamiokande detec- 
tor, at 2.5° relative to the proton beam direction. The on-axis 
detector (INGRID) consists of 7 +7 identical modules, arranged 
to form a "cross" configuration, and two "diagonal" modules 
positioned off the cross axes. The off-axis detector includes a 
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magnet, previously used in the UA1 and NOMAD experiments, 
operated with a magnetic field of up to 0.2 T; a Pi-Zero detector 
(POD); a tracking detector that includes time projection cham- 
bers (TPCs) and fine grained scintillator detectors (FGDs); an 
electromagnetic calorimeter (ECAL); and a side muon range 
detector (SMRD). 

The ND280 detector extensively uses scintillator detectors 
and embedded wavelength-shifting (WLS) fibers, with light de- 
tection from the fibers by photosensors that must operate in a 
magnetic field and fit in limited space inside the magnet. 

After studying several candidate photosensors, a multi-pixel 
avalanche photodiode operating in the limited Geiger multipli- 
cation mode was selected as the photosensor. These novel de- 
vices are compact, well matched to the spectral emission of 
WLS fibers, and insensitive to magnetic fields. Detailed in- 
formation about such devices and basic principles of operation 
can be found in recent review papers (see for example [3] and 
references therein). 

The operational parameters required for these photosensors 
were similar for all ND280 subdetectors and can be summa- 
rized as follows: an active area diameter of ~1 mm 2 , photon 
detection efficiency for green light >20%, a gain of (0.5 - 1 .0) x 
10 6 , more than 400 pixels, and a single photoelectron dark rate 
< 1 MHz. The pulse width should be less than 100 ns to match 
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the spill structure of the J-PARC proton beam. For calibra- 
tion and control purposes it was very desirable to obtain well- 
separated single electron peaks in the amplitude spectra for dark 
noise and low light levels. 

After an R&D study period of three years by numerous 
groups, the Hamamatsu Multi-Pixel Photon Counter (MPPC) 
was chosen as the photosensor for ND280. A description of 
this type of device and its basic parameters can be found in 
Ref. [4]. A customized 667-pixel MPPC with a sensitive area 
of 1.3x1.3 mm 2 was developed for T2K [5]. It is based on 
the Hamamatsu commercial device SI 0362-1 1-050C with 400 
pixels and lxl mm 2 sensitive area. The sensitive area was in- 
creased to provide better acceptance of light from 1 mm diam- 
eter Yll Kuraray fibers. In total, about 60,000 MPPCs were 
produced for T2K. The sensor is shown in Fig. 2. 
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Figure 2: Photographs of an MPPC with a sensitive area of 1.3 X 1.3 mm 2 : 
magnified face view (left) with 667 pixels in a 26 X 26 array (9 pixels in the 
corner are occupied by an electrode); the ceramic package of this MPPC (right). 

In this paper, we present the results of measurements and 
simulations of the main parameters of Hamamatsu MPPCs de- 
veloped for the T2K experiment, expanding upon the results 
given in [6]. We emphasize the operational parameters of these 
devices most critical for successful operation and calibration of 
the T2KND280 detectors: gain, dark rate, crosstalk, afterpulses 
and photon detection efficiency. This paper complements the 



results reported in [7], which focused on assessing the gross 
features of a large number of MPPCs. In this paper, dedicated 
setups were built to measure each process, which enabled more 
in-depth measurements than in [7] but in general these setups 
did not allow testing of a large number of MPPCs. 

2. MPPC response 

2.1. Operating principles 

A Multi-Pixel Photon Counter consists of an array of 
avalanche photo-diodes operating in Geiger mode. When op- 
erating in Geiger mode the diode is reverse-biased beyond 
the electrical breakdown voltage, which will be denoted Vbd 
throughout this document. Above Vbd, the electric field in 
the diode depletion region is sufficiently large for free carriers 
to produce additional carriers by impact ionization, resulting 
in a self-sustaining avalanche. In practice irreversible damage 
would eventually occur unless the avalanche is quenched. In 
MPPCs, quenching is achieved by using a large resistor in se- 
ries with the diode. The current produced by the avalanche cre- 
ates a voltage drop across the resistor (^quench), which stops the 
avalanche when the voltage across the diode reaches Vbd- The 
overvoltage, denoted AV, is the difference between the operat- 
ing voltage of the device and the breakdown voltage Vbd- The 
charge produced in an avalanche is hence the diode capacitance 
times AV. 

In Geiger mode, the amount of charge produced in an 
avalanche is independent of the number of charge carriers gen- 
erated within the depletion region. Hence, it is not possible 
to measure the light intensity by measuring the total charge 
produced in a single avalanche. MPPCs achieve photon count- 
ing capability by segmenting the detection area in an array of 
individual diode pixels. The amount of light hitting the de- 
vice is sampled by counting the number of pixels that pro- 
duce avalanches, which leads to a saturation effect when a large 
amount of light hits the sensor. However, the focus of this paper 
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is the MPPC response to low light levels, where the probabil- 
ity that multiple photons hit the same pixel at the same time is 
small. 

The T2K MPPC is an array of 26 by 26 pixels, each of which 
measures 50 x 50 //m 2 , on a common n+H — type silicon sub- 
strate [8]. Nine pixels in one corner have been replaced by a 
lead, reducing the total number of pixels to 667. The quench- 
ing resistors are polysilicon resistors. The Hamamatsu speci- 
fications sheet [4] states that the fill factor, i.e. the fraction of 
the device area that is active, is 61.5%. The breakdown voltage 
is about 70 V. When devices are purchased from Hamamatsu, 
rather than providing the breakdown voltage for each device, 
the voltage necessary to achieve a gain of 7.5 x 10 5 at 25°C is 
provided. 

2.2. Electrical properties 

The total resistance and capacitance of an MPPC were mea- 
sured using a picoammeter and capacitance-voltage (CV) an- 
alyzer, respectively. I — V and C - V plots are shown in 
Fig. 3. The MPPC capacitance was measured with a Keith- 
ley 590 CV analyzer. The capacitance drops rapidly with volt- 
age down to -20 V, which presumably corresponds to the full 
depletion of the device. The capacitance of the MPPC was 
found to follow a linear relationship when the supply voltage is 
less than -20 V: Cmppc = aV + b with a=0.0436+0.0003 pF/V 
and £=64.27+0.01 pF. At -70 V, the capacitance is then 
61.22+0.02 pF. The Hamamatsu specification document for 
T2K's MPPCs states that the terminal capacitance is 60 pF, 
which is consistent with 61.28 pF obtained at -70 V operating 
voltage. In the remainder of this paper, the minus sign will be 
omitted when discussing operating voltage. Using 60 pF total 
capacitance and neglecting parasitic capacitance yields a pixel 
capacitance of C p i x = 90.0 fF. 
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Figure 3: / - V and C - V plots for an MPPC. 

The current was measured with a Keithley 617 pro- 
grammable electrometer at 23°C. A linear fit for a for- 
ward bias voltage larger than 0.6 V yields a slope of 
^quench/(667 pixels)=225 O. From this we determine the av- 
erage quenching resistor value for this device to be /? qlle nch = 



150 kQ; for a set of thirty five sensors this parameter was dis- 
tributed in the range 148-154 kQ. 

2.3. Recovery time 

When an avalanche occurs in a pixel, the bias voltage across 
the diode drops down to the breakdown voltage. The diode volt- 
age recovers to the nominal operating voltage with a time con- 
stant that is nominally given by the product of the pixel capaci- 
tance and the quenching resistor. Using the values of /? quen ch and 
Cpi X ei reported in the previous section, the recovery time con- 
stant is t=13.4 ns. The overvoltage on the pixel at time t after 
the avalanche can then be written as: AV(f) = AV(0)(1 - e~'^ T ), 
where AV(0) is the nominal overvoltage. We will see in the fol- 
lowing section that the MPPC behavior is almost entirely driven 
by the overvoltage. Lower overvoltage implies a lower proba- 
bility of triggering an avalanche. It also implies a lower MPPC 
gain, hence an avalanche occurring while the pixel is recovering 
will yield a lower charge. 

The pixel voltage recovers to its nominal value by pumping 
charge from neighboring pixels and from the external electron- 
ics circuit. The capacitance of one pixel (90 fF) is small com- 
pared to the total capacitance of the MPPC (60 pF). Hence the 
voltage drop induced by the avalanche in one pixel on all the 
other pixels is very small. However, the neighboring pixels ef- 
fectively act as a bypass capacitor and the external circuit must 
eventually recharge the whole MPPC. The time constant intro- 
duced by the external circuit may be much longer than the pixel 
RC time constant and should be taken into account when in- 
vestigating the response of the MPPC to large light pulses, or 
when the repetition rate of avalanches is high. Since here we 
focus on characterizing the MPPC response to low light levels 
(<100 photoelectrons), the impact of the external electronics on 
the recovery time can be neglected. 

2.4. Photosensor gain 

The MPPC gain is defined as the charge produced in a sin- 
gle pixel avalanche, expressed in electron charge units. Single 
avalanches are typically created by a single carrier (unit charge) 
and can be triggered either by a photon or by thermal noise. 
Fig. 4 demonstrates excellent separation between the charges 
resulting from different number of photoelectrons. The gain 
is measured using a Multi-Channel Analyzer (MCA), by taking 
an amplitude spectrum and calculating the distance between the 
pedestal peak and the single photoelectron peak. Using other 
peaks provides consistent results. Conversion from the MCA 
output in digital counts to units of charge is achieved by cali- 
brating the electronics with a known input charge. The accuracy 
of the absolute gain measurement (i.e. the charge correspond- 
ing to a single avalanche) is affected mainly by the accuracy of 
the charge injection calibration. 

Fig. 5 shows the gain as a function of operating voltage for 
various temperatures. The curves are fit by linear functions ac- 
cording to gain, G = C p i x (V - VbdX where C p i x denotes the 
single pixel capacitance, V the operating voltage, and Vbd the 
breakdown voltage, i.e. the voltage for which extrapolation of 
the curves produces zero gain. The curves exhibit a slightly 
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Figure 4: A charge amplitude spectrum obtained using an LED source mea- 
sured with an MPPC (serial number TA9445) at an overvoltage of 1.5 V and 
temperature of 20° C. 

quadratic dependence, but a linear fit gives a reasonable esti- 
mate of Vbd and will be used throughout this paper. [We note 
that the voltage dependence of Cmppc reported in Section 2.2 
would cause the gain to have a quadratic dependence but this 
effect is smaller than the quadratic dependence we observe.] 
Since Vbd increases linearly by 52+4 mV per °C, the gain de- 
creases proportionately as the temperature increases at fixed op- 
erating voltage. However, the temperature variations within the 
T2K ND280 experiment are small enough that this effect can 
be calibrated out and does not require active compensation. 
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Figure 5: MPPC gain vs. supply (bias) voltage at different temperatures (sensor 
serial number TA8120). 

The overvoltage (AV) is calculated by subtracting the break- 
down voltage from the operating voltage. Fig. 6 shows the sin- 
gle avalanche charge as a function of AV. The fact that the 
curves lie on top of each other shows that the temperature de- 
pendence of the gain is dominated by the temperature depen- 
dence of Vbd- The slopes of the curves are consistent with the 
90 fF pixel capacitance estimated from the direct measurement, 



to within the equipment calibration accuracy. A detailed analy- 
sis of the temperature dependence of the capacitance measured 
at AV=1.3 V, shows a 0.1% increase per degree, which can be 
attributed to a change in the permittivity of the silicon [9]. 
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Figure 6: Single pixel charge (gain) of an MPPC (serial number TA8 1 20) as a 
function of the overvoltage AV at different temperatures. 

Fig. 4 shows that, unlike photomultiplier tubes, the MPPC 
gain fluctuations are significantly smaller than the charge from a 
single photoelectron avalanche. The gain fluctuations are, how- 
ever, not negligible. The spectrum presented in Fig. 4 can be 
fit by a series of Gaussian distributions, with the fi parameter 
for each Gaussian representing the mean charge in the peak and 
cr its width due to gain fluctuations and electronics noise. The 
gain fluctuation parameter cr(z') of the ;th peak is well described 
by the equation: 

°"0') 2 = O-ped + ' ■ ^Gain C 1 ) 

where cr pe d is the width of the pedestal, which is entirely due 
to the electronics noise, and cr Gain accounts for the gain fluctu- 
ations. Measurements of crGain show that it increases slightly 
with overvoltage. However, the achievable photoelectron res- 
olution is related to gain fluctuation relative to the measured 
gain, G, so in Fig. 7 we show the ratio 

G G w 

as a function of overvoltage, where <x pec i is the pedestal width 
and <x(l) is the width of the single avalanche peak. 

The 20° C data can be parameterized by the following func- 
tion: cr Gain /G = 0.064-AV~° 73 . The quality of the fit is good but 
we have no physical justification for this particular form. There 
appears to be a slight temperature variation, with the fluctua- 
tions being larger at higher temperatures. 

2.5. Dark noise 

Dark noise in Geiger-mode avalanche photodiodes is caused 
mainly by charge carriers generated thermally within the de- 
pletion region, which then enter the Geiger multiplication area 
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Figure 7: Relative gain fluctuation vs. overvoltage at various temperatures. The 
curve is a fit to the T=20°C data. 



and trigger avalanches. Any avalanche can, in turn, initiate sec- 
ondary avalanches through afterpulsing and crosstalk. Thus, 
the dark noise consists of single pixel avalanche pulses, along 
with larger amplitude pulses generated by optical crosstalk, af- 
terpulsing, and accidental pile-up from independent pixels. The 
last effect is negligibly small at dark rates below 1 MHz, as- 
suming a short integration time at the MPPC output. Optical 
crosstalk and afterpulsing are discussed in the next sections. 

Since most subsystems of our experiment acquire data as 
charge spectra within an integration gate associated with the 
beam crossing time, the relevant dark noise metric is the frac- 
tion of these gates populated by one or more dark pulses. The 
true rate of avalanches initiated by single charge carrier can be 
obtained from a Poisson distribution, using the following for- 
mula: 

/?DN = -ln(^)/A? (3) 

where no stands for the number of events with no counts, N 
for the total number of events, and Af for the gate time. The 
measurements presented here used 160 ns gates triggered at a 
constant rate of 20 kHz. 

Fig. 8 shows that the dark rate increases linearly with over- 
voltage in the range of 0.5-1.6 V. Above 1.6 V the points devi- 
ate upwards from the linear fit, which we attribute to an effect 
of afterpulsing. The temperature dependence is exponential and 
is shown in Fig. 9. The data for each sensor has been fit with a 
function of the form given in Eq. (4). 
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where T is absolute temperature. In this formula A represents 
the ratio of dark rate to overvoltage at T=298 K (25°C) (in 
kHz/V). Vo is the offset of breakdown voltage calculated from 
the dark rate with respect to that obtained from the gain mea- 
surements, and E the band gap energy. The fit range was re- 
stricted to AV < 1.6 V and R DN < 5 MHz, in order to avoid the 
effect of afterpulsing and rate limitations of the equipment. 



The parametrization given in Eq. (4) has been obtained under 
following assumptions: 

1 . A non-degenerate semiconductor model was used. 

2. Thermally generated charge carriers are a result of trap- 
assisted (i.e. involving an R-G center 1 ) generation pro- 
cesses. 

3. Given high reverse bias, the device operates in the so 
called 'R-G depletion region' steady state, i.e. no free 
charge carriers are available within the depletion region. 

4. The trap energy level is close to the middle of the silicon's 
bandgap. 

5. Only processes occuring within the volume of the de- 
pletion region are taken into account. Surface genera- 
tion/recombination is neglected. 

Using such model, one can easily explain significant sensor-to- 
sensor variations of the dark rate by: a) differences in the con- 
centrations of traps (R-G centers) and b) differences in dopant 
concentrations, hence different junction volumes. Mean value 
of the observed bandgap energy for the five measured sensors 
(Table 1) is 1.127+0.0099 eV, which is within the range of val- 
ues widely reported for silicon. Furthermore, reasonable x 2 / v 
values and an average /?-value of 33.4% do not provide enough 
evidence to reject the parametrization given by Eq. (4) at a sta- 
tistically significant level, which is why we assumed that it can 
be used to approximate data from our measurements. 




2.5 
AV (V) 

Figure 8: Dark rate vs. overvoltage AV at different temperatures (sensor num- 
ber TA8120). A single fit to Eq. (4) has been used to fit all the data points. Solid 
lines show results within the fit range (AV < 1.6 V and Sdn 2 5 MHz) while 
dashed lines represent extrapolations. 

The dark noise rate varies significantly between MPPCs as 
reported in [7]. A 20% variation in the dark noise rate was 
found at 20°C and AV = 1 V when calculating the variation 
as the root mean square (RMS) over the mean for the 17,686 
tested MPPCs. 
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Figure 9: Dark rate vs. temperature at different overvoltages AV (sensor num- 
ber TA8120). A single fit to Eq. (4) has been used to fit all the data points, with 
the fit range restricted to AV < 1.6 V and R D N < 5 MHz. 



Table 1: Fit parameters for the dependence of the dark rate on AV and tem- 
perature. Eq. (4) was used with the fit range restricted to AV < 1.6 V and 
Rdn < 5 MHz. A is the dark rate to overvoltage ratio at T=298 K (25°C). Vo 
is the difference between the breakdown voltage calculated from gain and dark 
noise. E is the band gap energy. 



Sensor no. 


A (kHz/V) 


V (mV) 


£(eV) 


2 i 


TA8744 


934+5.0 


73±5.2 


1.117±0.0024 


1.10 


TA8160 


562+3.7 


98+6.1 


1.139+0.0030 


0.99 


TA8120 


564±3.7 


93+6.4 


1.135+0.0035 


1.15 


TA8092 


622±3.8 


92+5.9 


1.126+0.0034 


1.07 


TA9314 


789+3.8 


72+4.7 


1.118+0.0022 


0.71 


Mean 






1.127 ±0.0099 





2.6. Afterpulsing 

2.6.1. Correlated noise 

Correlated noise is a general label for avalanches that are 
triggered by other avalanches. There are two known types of 
correlated noise: crosstalk and afterpulsing, both of which will 
be described and characterized in details in the next two sub- 
sections. In general, whenever an avalanche occurs there is a 
chance that it triggers one or more additional avalanches, either 
in neighboring pixels or in the same pixel at a later time. As 
mentioned earlier, the dark noise rate was estimated by count- 
ing the number of time no avalanches were detected within a 
gate. Indeed, the average number of avalanches detected within 
the gate is not a good estimator of the dark noise rate because 
some avalanches may have occurred due to correlated noise. 
Hence, in the presence of correlated noise, the measured av- 
erage number of avalanches will exceed the expectation from 
Poisson statistics. Conversely, the measured number of events 
having one avalanche within the gate will be lower than the 
expectation. This fact can be used to get an estimate of the 
correlated noise. 

The data used for measuring dark noise presented in the pre- 
vious section can also be used to get an estimate of the cor- 



related noise. From the dark noise rate measurement one can 
predict the fraction of events that should have one avalanche 
in the absence of correlated noise. The correlated noise proba- 
bility fcN is the probability that one avalanche triggers at least 
one additional avalanche. The presence of this correlated noise 
term modifies the calculation of the number of events with one 
avalanche within the gate, N\, as follows, 



= e- R ™ A 'R DN At(l-P CN ) 



(5) 



where iV tot is the total number of events recorded, /?dn is the 
dark noise rate, and At is the gate width. The correlated noise 
probability estimated by solving this equation for Pcn is shown 
in Fig. 10 at different overvoltages and temperatures. 
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Figure 10: Combined crosstalk and afterpulse probabilities vs. overvoltage at 
several temperatures. 

The temperature dependence is strikingly small. The data 
can be fitted by a quadratic function: Pcn = kAV 2 with 
k=9.4+0.1 %. The quadratic fit is good until AV > 1.6 V, which 
is also approximately the overvoltage at which the measured 
dark noise rate no longer behaves linearly. As explained earlier, 
at sufficiently large overvoltage the method used to estimate 
dark noise becomes compromised by afterpulse avalanches that 
stem from dark noise avalanches prior to the integration gate. 
Hence, it is likely that the failure of the fit results from the dark 
noise rate being incorrectly inferred when AV is greater than 
about 1.6 V. 

2.6.2. Measuring afterpulsing 

Afterpulsing is understood as being caused by the trapping 
of charge carriers created during an avalanche. The trapped 
carriers eventually get released and trigger an avalanche within 
the same pixel as the original avalanche, but delayed in time. 
Afterpulsing may be partially suppressed by the fact that the 
pixel voltage recovers in about 45 ns (a 13.4 ns time constant) 
as described in Section 2.3. If a carrier is released while the 
pixel voltage has not reached the nominal voltage, then the 
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charge produced in the avalanche will be lower than for nomi- 
nal avalanches. For a self-consistent description of the data, it 
is best to factorize recovery and afterpulsing phenomena, i.e. 
to measure the number of afterpulse avalanches per original 
avalanche regardless of the pixel voltage at the time of the after- 
pulse avalanche. Because there may be different types of traps 
in the silicon, there is no reason to assume that afterpulsing 
should follow a single exponential decay. In fact previous mea- 
surements on a similar MPPC [10] have shown that afterpulsing 
exhibits two time constants. 

Two methods were used to measure afterpulsing. The meth- 
ods complement each other since they effectively probe dif- 
ferent afterpulsing time constants. Both rely on the fact that 
afterpulse avalanches are correlated in time with their parent 
avalanche. 

The first method is based on the analysis of waveforms de- 
scribed in [10]. The waveforms were fit with a superposition 
of single avalanche response functions that allow separating 
pulses occurring within a few nanoseconds of one another. The 
probability of an avalanche occurring at time t after another 
avalanche can be parameterized as: 



P(t) 



1 



f 

Jo 



P AP (x)dx -P DN (t) + 1- P DN (x)dx -P AP (t) 

L Jo 

(6) 

where Pap and Pdn are the afterpulsing and dark noise prob- 
abilities. [We note that the formula in Ref. [10] is incorrect 
and should be replaced by this one.] Afterpulsing can be pa- 
rameterized using two parameters: «ap, the number of after- 
pulse avalanches created per original avalanche, and t, the time 
constant of the exponential distribution governing the afterpulse 
generation. A drawback of this method is that the likelihood of 
uncorrelated dark signals (with a rate of 500-1000 MHz) in the 
waveform limits the sensitivity to afterpulsing time constants of 
less than about 100 ns. However, this method is well suited for 
measuring small time constants (< 50 ns) as the pulse finding 
techniques allow detecting pulses separated by a few nanosec- 
onds. 

The second method is based on counting the number of 
avalanches with a scaler after introducing a controlled deadtime 
following each detected avalanche. The width of the analog 
pulse resulting from the convolution of the MPPC and ampli- 
fier response was such that the minimum deadtime that could be 
set was 26 ns. This minimum gate width means that this method 
is sensitive only to afterpulsing time constants of greater than 
about 50 ns. However, in contrast to first method, the count- 
ing technique overcomes the dark noise background limitation 
when measuring long time constants by taking very high statis- 
tics data. The deadtime dependent rate can be fit by a func- 
tion that includes the contribution of dark noise and afterpuls- 
ing. In the absence of afterpulsing, the measured rate R(At) 
for a given deadtime At is Pdn/(1 + PdnA?)- Afterpulsing pro- 
duces avalanches that will increase the rate as long as they oc- 
cur after the deadtime. To first order (i.e. assuming that one 
avalanche creates at most one additional detectable afterpulse 
avalanche and ignoring afterpulse avalanches created by previ- 
ous afterpulse avalanches) the measured rate can then be calcu- 



lated from: 

R(At) = R/(l +RAt) (7) 
with (assuming two afterpulsing time constants) 

R = Wd " n APO e- A " T " - n APl e- A " Tl ) (8) 

where /?dn is the dark noise rate, «APi = 0, 1) is the average 
number of afterpulse avalanches per original avalanche, and t, 
(; = 0, 1) the afterpulsing decay time constant. 




Figure 1 1 : Afterpulse parameters (exponential time constant and probability of 
having at least one afterpulse) vs. overvoltage for four different MPPCs (with 
plot symbols: square, circle, triangle, diamond) at 25°C. 

Fig. 11 shows the afterpulse parameters for four different 
MPPCs measured at 25°C using the waveform analysis tech- 
nique. The probability is calculated from the number of after- 
pulse avalanches per original avalanche distributed as 1 - e~" AP 
and so is the probability that an avalanche generates at least 
one afterpulse avalanche. The exponential time constants were 
found to be r s hoi-t = 17.6 ± 2.1 ns and ri ong = 71.4 + 8.3 ns. 
The probabilities of "long" and "short" afterpulses are almost 
equal. The total probability of afterpulses is about 0.16 per 
initial avalanche at AV=1.4 V. The number of afterpulses per 
avalanche as a function of overvoltage can be fit by a sim- 
ple quadratic function: n A p(AV) = K ■ AV 2 , with A'short = 
0.0400±0.001(sfaf)+0.005(svs)y- 2 and Jr r long = 0.0402(sfaf) ± 
0.001 + 0.005(sy.s)V~ 2 . The dominant systematic error arises 
from the inability to detect pulses less than 2 ns after the first 
avalanche. 

Fig. 12 shows the rate measured as a function of deadtime 
from 26 ns to 1 //s at an overvoltage of 1.4 V at 25°C. The fit to 
the data is excellent with an average^ 2 of 74.8 for 95 degrees of 
freedom. The fit parameters for the data from the MPPC shown 
in Fig. 12 yield t = 57+5 ns, n A p = 0.107 + 0.005,^ = 287± 
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Figure 12: Rate as a function of deadtime for an MPPC biased at 1.4 V over- 
voltage. The black curve shows a fit function including dark noise and two 
afterpulsing time constants. The dashed curve shows the estimated contribu- 
tion of dark noise, i.e. the fit to Eq. (8) with afterpulsing probabilities turned 
off. The upper right inset shows the same data after subtracting the estimated 
dark noise. 



49 ns, and «api = 0.043 + 0.006. Repeating this test over 35 
different MPPCs yield the following averages: to = 52(8) ns, 
«apo = 0.105(0.009), T\ = 315(84) ns and n A?l = 0.066(0.01), 
with the standard deviations indicated in parentheses. 

Since the two measurement methods are sensitive to differ- 
ent afterpulsing time constant ranges it is not surprising they 
yield different results. It is possible to reconcile both meth- 
ods by fitting the variable deadtime data for all MPPCs with 
three different afterpulsing time constants, two of them be- 
ing fixed: tq-\1 ns and ti=70 ns. The third time constant 
is a free parameter in the fit. Excellent fits are again ob- 
tained with an average^ 2 of 72.52 for 95 degrees of freedom, 
which is slightly better than the fit with just two time con- 
stants. The parameters averaged over the 35 MPPCs are «apo = 
0.058(0.03), hapi = 0.090(0.008), n AP3 = 0.056(0.009), and 
T3 = 373(55) ns, with the standard deviation among MPPCs 
in parentheses. The expectation from the waveform analysis at 
1.4 V is «AP s hon = 0.078 and «APi on o = 0.082. The introduction 
of this third (373 ns) time constant into the fitting function used 
for the waveform analysis at the level suggested by the vari- 
able deadtime analysis does not worsen the agreement with the 
data significantly and so is an acceptable additional parameter 
in a range not accessible to the method. Both analyses are also 
qualitatively consistent with the simple correlated noise analy- 
sis presented in the previous section, which predicts a total con- 
tribution of 0. 184 for crosstalk and afterpulsing at AV = 1 .4 V. 

The temperature dependence of afterpulsing was measured 
with the waveform technique for a couple of MPPCs at con- 
stant AV within a range of 13-25°C as shown in Fig 13. The 
amplitude of the "long" component of the afterpulsing rate is in- 
sensitive to temperature within measurement accuracy. On the 
other hand, the amplitude of the "short" component decreases 
as the temperature increases with a coefficient of 2.0-2.5% per 
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Figure 13: Afterpulse parameters (exponential time constant and probability of 
having at least one afterpulse) vs. temperature. 



°C. The short and long time constants decrease with increasing 
temperature from 21 ns and 90 ns at 13°C to 17 ns and 70 ns 
at 25°C respectively. In the MPPC simulation code, the tem- 
perature dependence of the 17 ns and 70 ns time constants is 
implemented, while the 370 ns time constant will be assumed 
to be constant since no temperature dependent data are available 
to quantify the variation. 

2. 7. Optical crosstalk 

2.7.1. Crosstalk measurement 

Optical crosstalk is believed to occur when optical pho- 
tons produced in an avalanche propagate to neighboring pix- 
els where they produce photoelectrons [11]. The result is that 
two or more pixels can be fired almost simultaneously (i.e. on 
a timescale of «1 ns). The photon emission probability has 
been estimated to be 10~ 5 photons per carrier crossing the junc- 
tion [12], the absorption length for photons that are most effec- 
tive in propagating the avalanches is typically 1 mm. Although 
the total crosstalk fraction is small, it is expected to vary with 
overvoltage and a detailed study is necessary to fully charac- 
terize local variations of the crosstalk phenomenon. Measure- 
ments of crosstalk variations within the pixel array as a function 
of voltage were performed using optical microscopy and wave- 
form analysis. 

2.7.2. Optical microscopy 

Crosstalk probabilities were measured using the apparatus 
shown schematically in Fig. 14. A nanoLED [13] light source 
system was used to produce a pulse width of 1 ns FWHM from 
an integrated 463 nm LED. This was coupled to an optical fiber 
that terminated in a microscope lens such that the light beam 
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Figure 14: Schematic of the optical microscopy apparatus used to measure 
MPPC crosstalk within a pixel. 



is focused onto the MPPC face. The MPPC was mounted on 
an X-Y stage so that the beam spot could be translated across 
the MPPC pixel array with one micron position resolution. The 
MPPC signal was digitized with a 1 GHz sampling rate during a 
1 fis period using a Tektronix TDS 380 oscilloscope. The light 
pulse intensity was measured to be between 10-20 photons at 
the exit point of the fiber. 

It was assumed that the amplitude of the avalanche signal ob- 
served is not dependent on the number of photons injected into 
a pixel if the photons all originate from the same LED pulse. 
Therefore each LED flash creates a 1 p.e. signal as long as the 
beam spot is well-centered within a pixel. The trigger efficiency 
was measured as the beam was scanned across several pixels 
to estimate the profile of the photon beam. The profile mea- 
sured over 150 /vm is shown in Fig. 15 and is consistent with a 
Gaussian beam spread of 5 fim. The sensitive area is in good 
agreement with the value of 61.5% specified in the Hamamatsu 
catalog. 
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Figure 15: MPPC detection efficiency across several MPPC pixels. The dashed 
line is the expected profile from the convolution of a Gaussian beam spot with 
<x = 5 /im and square 50 /im pixels. 



2. 7.3. Crosstalk study using waveform analysis 

Waveforms were recorded for nine beam positions inside 
each of three pixels chosen for their specific position within the 
array, namely: a corner pixel, a side pixel and a pixel inside the 
MPPC array away from the edges. These pixels are surrounded 



by 3, 5 and 8 pixels respectively, each of which may generate 
crosstalk signals. Crosstalk probabilities were calculated for in- 
dividual photoelectron pulses selected to be within 8 ns after the 
LED trigger pulse. The total crosstalk probability is given by: 



Pot = 1 - 



jV(lpe) 
Mot ' 



(9) 



where N(lpe) is the number of single 1 p.e. pulses and N tot 
is the number of all LED pulses. The total crosstalk signal is 
defined as the observation of >2 p.e. pulses within the 8 ns time 
window, while individual crosstalk probabilities were extracted 
by selecting pulse heights corresponding to 2 p.e., 3 p.e. and 
4 p.e. Data was taken for overvoltage AV = 1.335 V and T = 
24°C; results for all three pixels are presented in Fig. 16. For all 



Corner Pixel (Trigger pulse) 



0.045 0.077 0.087 



0.037 0.055 0.052 



0.085 0.091 0.108 



0.122 0.116 



0.112 0.111 0.122 



0.126 0.109 0.133 



Figure 16: Measured crosstalk probabilities for nine beam positions inside an 
MPPC pixel at AV = 1.335 V. (Left) a corner pixel with the corner located on 
the bottom left; (center) a side pixel with the side boundary to the left of the 
pixel and (right) a pixel inside the MPPC array. The thick black line denotes 
the limit of the pixel array. 

three pixels the crosstalk measured shows a clear dependence 
with position of the beam spot, suggesting that the crosstalk 
probability is correlated with where the photon is absorbed in 
the pixel. A similar analysis was applied to MPPC dark count 
data in a time window 500 ns before the LED triggers. The 
crosstalk was measured to be 9 + 1% and no correlation with 
beam spot location was found. 

Fig. 17 presents measurements of crosstalk probabilities as 
a function of overvoltage at T — 24°C for the same three pix- 
els. Based on the position variation results, a correction factor 
is applied to correct the crosstalk probability to a probability 
averaged over the entire pixel. All probabilities were found to 
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Figure 17: Crosstalk value vs. overvoltage for three pixels shown in Fig. 16. 
P a is crosstalk probability, P ct =i is the probability of only one pixel fired in 
addition to the initial pixel, etc. 

agree with a AV 2 dependence except for the corner pixels where 
the total probability plateaus at high overvoltage (above 1.3 V). 
This plateau is due to some peculiar behaviors of the corner pix- 
els, which cannot be explained by geometrical considerations. 
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Variations of the total crosstalk probability between pixels is in 
good agreement with the hypothesis of a point source genera- 
tion of optical photons in the pixel. The result of the crosstalk 
simulation is shown in Fig. 18 and agrees well with the data for 
a nearest neighbor crosstalk hypothesis. This model is included 
in the simulation described in Section 3. 
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Figure 18: Total crosstalk probabilities for all three pixels (corner, side and 
in-array). Data is shown as a solid circle, a solid black line indicates the best 
quadratic fit from Fig. 17. Simulated probabilities are shown as open circles 
and a dashed line. 



2.8. Photon detection efficiency 

The photon detection efficiency (PDE) of a multi-pixel 
avalanche photodiode operated in a limited Geiger mode is a 
product of three factors: 



PDE = QE(A) ■ eceiger ■ s pixe i, 



(10) 



where QE(A) is the wavelength dependent quantum efficiency, 
fiGeieer is the probability to initiate the Geiger discharge by a 
carrier, and geometric acceptance e p i xe i is the fraction of the 
total photodiode area occupied by the photosensitive area of the 
pixels. 

For an MPPC, quantum efficiency can be defined as the prob- 
ability for an incident photon to generate an electron-hole pair 
in a region in which carriers can produce an avalanche. The 
layer structure in an MPPC is optimized to have the highest 
probability for a visible photon to be absorbed in the depletion 
layer. For comparison, an APD with a similar layer structure 
to that of the MPPC developed by Hamamatsu for the CERN 
CMS experiment [14] has a measured quantum efficiency of 
more than 80% at 500 nm, so a similar value may be expected 
for the MPPC. 

Overvoltage affects just one parameter in the expression, 
namely fiGeiger. The breakdown probability depends on the 
impact ionization coefficients (for electrons and holes), which 
are strong functions of electric field. Simulation and measure- 
ments [15] show that SGeiger behaves as the exponentially sat- 
urating function s ma x(l - e~ kAV ) if breakdown is triggered by 
electrons. Breakdown initiated by holes leads to a linear depen- 
dence on AV. 



The geometrical factor e p i xe i is solely determined by the 
MPPC topology. Our measurements indicate g p j xe i=64%, which 
is consistent with the Hamamatsu specification of 62% for sen- 
sors with 50 fim pixels. 

2.8.1. PDE measurement-pulsed LED method 

For the PDE measurements we used an approach discussed 
in [16]. The PDE is measured using pulsed LED light with a 
narrow emission spectrum. The number of photons per LED 
flash is collimated to be within the MPPC sensitive area and 
reduced to an intensity that can fire only 2-5 pixels on average. 
The number of photons per LED pulse N y (A) can be measured 
using a calibrated photodetector, i.e. one with known spectral 
and single electron responses. 

The PDE can be calculated from the recorded MPPC pulse 
height distribution (see Fig. 4) by assuming a Poisson distri- 
bution of the number of photons in an LED pulse. The mean 
value jVp e of the number of photons recorded per LED pulse 
can be determined from the probability P(0) of the pedestal 
(0 p.e.) events by N pe = -lnP(0). N pe calculated this way is 
independent of afterpulsing and crosstalk. Knowing the num- 
ber of photons incident on the MPPC, N 7 (A), one can calculate 
PDE(A) = N P JN Y (A). 

The dependence of PDE on bias voltage was measured using 
a fast green emitting LED operating in a pulsed mode. The 
emission spectrum of this LED was measured to be very close 
to that of a Y 1 1 WLS fiber. The peak value is centered around 
515 nm, and it has FWHM of 40 nm. 

The MPPC was illuminated with LED flashes through a 
0.5 mm collimator (the distance between the LED and the colli- 
mator was 20 cm, and with about 1 mm between the collimator 
and the MPPC). A neutral density filter reduces the light inten- 
sity on the MPPC face to the level of 10-15 photons. The signal 
from the MPPC was amplified with a fast transimpedance am- 
plifier and digitized with a Picoscope 5203 digital oscilloscope 
(250 MHz bandwidth, 1 GHz sampling rate). The signal inte- 
gration time was 150 ns. A schematic diagram of the setup is 
shown in Fig. 19. 

Picoscope 5203 
gate 




200 mm i 1 mm 



Figure 19: Schematic diagram of the setup for the PDE measurements. 

The PDE measurements were done in a temperature stabi- 
lized dark box (AT<0. 1 °C). The number of photons in the LED 
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flash was measured using a calibrated Hamamatsu photomulti- 
plier R7899 (QE=15.7% at 515 nm). The photoelectron collec- 
tion efficiency for the 5-mm diameter central part of the photo- 
cathode is more than 95%, the PMT excess noise factor is 1.15. 
The LED amplitude spectrum measured for one of the tested 
MPPCs is shown in Fig. 4 at AV=1.5 V and 20°C. 

The average number of photoelectrons in the LED signal was 
calculated by counting the number of pedestal events as dis- 
cussed in Section 2.5. To correct for dark pulses that occurred 
randomly inside the 150 ns integration gate, the dark rate was 
measured during the same gate width but 300 ns earlier relative 
to the LED pulse. The stability of the LED pulse intensity was 
monitored and found to be better than +3% during the measure- 
ments. 

MPPC PDE as a function of overvoltage is shown in Fig. 20 
for three temperatures. The PDE depends almost linearly on 
AV within the AV range of 1.0-1.6 V with slope of 1.5% per 
100 mV. For a fixed overvoltage there is no observable depen- 
dence on temperature. 
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Figure 21: PDE (at 515 nm) vs. overvoltage for four MPPCs at 20°C. 
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Figure 20: Photon detection efficiency of a MPPC (serial number TA9445) for 
green light (~515 nm) as a function of overvoltage at three temperatures. 

Fig. 21 shows the measured PDE for four additional MPPCs 
at 20°C. All show essentially the same performance with the 
PDE in the range 29-32% for green light at a typical operating 
overvoltage of 1.4 V. The measurement accuracy of the PDE 
is estimated to be about 10%. The largest contribution to this 
uncertainty is the normalization error, which is dominated by 
the error in the PMT spectral calibration (5-7%) followed by the 
uncertainty in the p.e. collection efficiency in the PMT (5%). 

2.8.2. PDE measurement-optical power meter method 

The PDE has also been measured independently using a 
473 nm LED pulser developed as a calibration source for the 
ANTAPvES experiment [17] and with a 463 nm NanoLED 
pulser. The experimental setup is similar to that shown in 
Fig. 19. The MPPC signal was amplified with a gain of 40 
and then sampled by a LeCroy Wave Runner 6100 oscilloscope 
(1 GHz bandwidth, 10 GSample per sec) within a 200 ns 



gate. The temperature was held stable to 0.1 °C by means of 
a thermally-coupled metallic plate. 

Calibration of the number of photons incident on the MPPC 
was performed using a Newport 1835-C Optical power meter 
with an 818 series 1 cm 2 photodiode sensor. The power meter 
converts an optical signal of specific wavelength into the opti- 
cal power equivalent. The number of photons per flash can be 
calculated as follows, 



P w -A 
iVv = — — - A, 



'Hz 



he 



(11) 



where Pw is the measured optical power (in watts) at wave- 
length A (463 or 473 nm), F Hz is the LED pulse rate (13 kHz), h 
is Planck's constant and A is an acceptance factor. Acceptance 
A is the ratio of intensity of LED light incident on the MPPC 
sensitive area to the intensity incident on the power meter sen- 
sor. The value of A is evaluated by performing a position scan 
of the LED light intensity profile. 

The power meter calibration accuracy and the estimated ac- 
ceptance factor are the dominant contributors to the uncertain- 
ties in the PDE measurement. As discussed in Section 2.8.1, 
the number of detected photoelectrons was obtained from the 
number of pedestal events in signal and assuming a Poisson dis- 
tribution. The PDE values are free of afterpulsing and crosstalk 
contributions and were corrected for dark rate. The PDE was 
found to be 3 1 % for the 463 nm LED and 29% for the 473 nm 
LED at AV=1.3 V, which is in good agreement with the MPPC 
spectral sensitivity shown in Fig. 22 and discussed in the next 
section. 

2.8.3. Spectral sensitivity 

A spectrophotometer calibrated with a PIN-diode [16] was 
used to measure the spectral sensitivity of the MPPC. The spec- 
trophotometer light intensity was reduced until the maximum 
MPPC current was only -30% greater than the dark current 
to avoid nonlinearity effects caused by the limited number of 
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pixels. Comparing the MPPC current with the calibrated PIN- 
diode photocurrent we obtain the relative spectral sensitivity. 

To achieve an absolute scale, the measured relative spectral 
response is scaled to the reference PDE points obtained with 
LED light at three wavelengths: 410, 460 and 515 nm mea- 
sured at 1.2 V overvoltage. The scaling factor at other overvolt- 
age values was found to be constant at these wavelengths up to 
about 1 .4 V — the PDE spectral sensitivity shape is appreciably 
different above this, as was noted in [15]. The MPPC PDE de- 
pendence on the wavelength of the detected light along with the 
emission spectrum of the WLS Yl 1 Kuraray fiber are shown in 
Fig. 22. The MPPC peak sensitivity is in the blue light region, 
around 450 nm. 

Since the spectrum of light incident on the MPPC in the 
ND280 detectors is determined by the Yl 1 fiber emission spec- 
trum and the wavelength-dependent attenuation in the fiber, a 
PDE measurement was performed by exciting a Yl 1 fiber with 
a 405 nm LED [18]. The blue light source was arranged so that 
only re-emitted green light could reach the photosensor after 
propagating through 40 cm of the fiber. The fiber was coupled 
directly to the MPPC with the same design of optical connector 
used in the ND280 ECAL and POD subsystems. At AV=1.3 V, 
the PDE was measured to be 21%, which is significantly lower 
than the 28% measured at the same overvoltage with light in- 
cident directly onto the MPPC. The lower value may be due to 
light loss at the interface between the coupler and the Yl 1 fiber. 

Fig. 22 also shows the MPPC spectral sensitivity measured 
by Hamamatsu for a commercial MPPC SI 0362-1 1-050 device 
at 25°C. These data, taken from the Hamamatsu catalog, are not 
corrected for crosstalk and afterpulsing. The method Hama- 
matsu used is basically the photocurrent method described 
above but with a monochromator to select the incident light 
wavelength. The number of incident photons is derived from 
a calibrated photodiode response and the number of detected 
photoelectrons is obtained by dividing the MPPC current by its 
gain and the charge on an electron and assuming a Poisson dis- 
tribution of the number of photons per single flash. We have 
corrected the Hamamatsu result by scaling down the PDE val- 
ues by 0.663. This scaling factor was chosen to fit the sensitiv- 
ity curve at the points measured with LEDs; the renormalized 
Hamamatsu spectral plot shape is consistent with our results 
within measurement accuracy. 

3. MPPC simulation 

A Monte Carlo simulation of the MPPC, written in C++, is 
now at a mature stage of development. The simulation can be 
split into two main components - a set of models defining device 
behavior, and a procedural framework to initialize the model 
using input parameters, control the simulation and output the 
results. The framework will be discussed briefly first. 

3.1. Simulation framework 

The simulation is based on a list of potential trig- 
gers (incident photons, thermally generated carriers and 
crosstalk/afterpulses), which are processed in time order. The 
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Figure 22: MPPC PDE as a function of wavelength at AV=1.2 and 1.5 V at 
25°C. Also shown is the spectral plot from the Hamamatsu catalogue, which 
uses data not corrected for crosstalk and afterpulsing (blue line); the dashed line 
is the Hamamatsu plot scaled-down using knowledge of the correlated noise 
contribution from our measurements. The green curve shows the Yl 1(150) 
Kuraray fiber emission spectrum (arbitrary units) for a fiber length of 150 cm 
(from Kuraray spec). LED and Yl 1 fiber points were measured at AV=1.3 V. 



only state variables of the MPPC are the voltages across each 
pixel; the evolution of these voltages between triggers is han- 
dled by a recovery model. On initialization, a list of incident 
photons is given to the simulation as input, and thermal noise is 
generated at the appropriate rate DCR(V nom ) for a nominal op- 
erating bias voltage and temperature. These two sources form 
the initial list of potential triggers. 

Each potential trigger is then processed in the following 
steps: 

1 . The voltages on all pixels are updated from their state af- 
ter the previous trigger, using the recovery model and the 
elapsed time since the last processed trigger. 

2. It is determined whether the pixel fires. The probability is 
equal to PDE( V pix ) for photons and DCR( V pix )/DCR( V nom ) 
for dark noise, to account for the lower DCR for a pixel 
with depleted voltage, relative to the nominal DCR used 
to generate the noise triggers. 

3. If the pixel fires, a trigger is added to the list of output 
signals and its voltage is set to zero; the charge of the 
generated avalanche depends on the voltage of the fired 
pixel and it is smeared by a Gaussian resolution function. 
The afterpulse/crosstalk models determine whether further 
noise is generated, and, if applicable, the additional noise 
is inserted into the list of potential hits, in correct time or- 
der. 

The reinsertion of correlated noise resulting from an initial 
trigger allows higher-order noise cascades to be dealt with in a 
simple and natural way. The final output is a list of avalanches 
with times and charges, which can then be processed by code 
appropriate to a specific readout circuit. 
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3.2. Physics models 

The simulation relies on accurate models for the various ef- 
fects present in the sensor. The characterization measurements 
described above have been used to determine appropriate mod- 
els to use, and to tune model parameters. 

The dark rate is parameterized as a linear function of bias 
voltage-the parameters for this function must be measured sep- 
arately for each sensor since large variations between devices 
are observed. The PDE is modeled with a quadratic fit to the 
data in Fig. 20; variation with wavelength is not included. The 
mean number of short- and long-lived trapped carriers for af- 
terpulsing, and the lifetime of the trapped states, are taken from 
the results of the waveform analysis method in Section 2.6. The 
crosstalk model is based on the data and the model described 
in Section 2.7. The data shown in Fig. 18 are well-described 
by a simple nearest-neighbor model that assumes crosstalk oc- 
curs only in the four nearest-neighbor pixels to the primary 
pixel. Crosstalk pulses are generated according to the proba- 
bility measured from dark noise. The location of the crosstalk 
pulse is then chosen randomly among the four neighboring pix- 
els. The pulse is discarded if it falls outside the MPPC active 
area. 

The recovery model used is specific to the ND280 Trip-t- 
based electronics (TFB board [19]), it assumes recharging of 
the fired MPPC pixels from capacitances elsewhere in the read- 
out electronics for each channel. Recovery does not signifi- 
cantly affect response at low light levels however, so it will not 
be discussed in further detail. 

3.3. Comparison with data 

The simulation output has been compared to data taken us- 
ing the ND280 Trip-t electronics and a fast-pulsed LED, with 
a gate length of 540 ns and the photosensor at a temperature 
of 22°C. An adjustable lens was used to alter the intensity of 
light incident on the sensor. All the parameters used for the 
simulation were taken from the characterization measurements, 
but some tuning was required to reflect sensor-specific parame- 
ters, electronics effects and light-level uncertainties. The linear 
fit parameters for the sensor dark noise curve were measured 
and used in the simulation. Since an absolute calibration of the 
incident light level was not available, the mean incident pho- 
ton number was calculated for 1.33 V overvoltage using the 
method described in Section 2.8.1. The absolute PDE in the 
simulation is therefore not tested by this comparison, but errors 
in the parametrization of PDE as a function of voltage will be 
evident. Finally, the spread in total event charge due to electron- 
ics noise, and the spread in avalanche gains, were determined 
from the measured peak widths at a low light level and 1.33 V 
overvoltage, and added to the simulation. 

Histograms of integrated output charge are shown in Fig. 23, 
for data and simulation. Very good agreement is seen between 
the data and MC for a range of light levels and overvoltages. 
Some small discrepancies between data and MC are seen in 
the integer-binned histograms; however these histograms de- 
pend on the exact peak positions, which must be determined in 
the data by fitting. They also depend sensitively on the exact 




Figure 23: Comparison of data to Monte Carlo at low light level for a range of 
overvoltages. The photon numbers shown are the number incident on the face 
of the MPPC. The histograms on the right show the same data as on the left, 
but with a bin width equal to the fitted peak separation in the data. 

shapes of the peaks, since for large peak widths, each integer 
bin contains some events which have migrated from neighbor- 
ing peaks. No significant systematic difference is observed be- 
tween data and MC. 

3.4. Energy resolution 

In most cases, the energy resolution of scintillator detectors is 
dominated by the photon counting statistics when the number 
of photoelectrons is low (less than about 100). However, the 
photosensor and electronics can impact the energy resolution 
in two ways: 1) constant noise background due to dark noise 
and electronics noise, 2) fluctuation in the charge detected per 
photoelectron. The energy resolution can be calculated fairly 
accurately in the case where the MPPC charge is integrated over 
a time window At and ignoring the MPPC saturation effect. The 
standard deviation of the number of avalanches can be written 
as: 

°n*, = ^Av + A^AvO-Q + °i + ^dnA? (12) 

where Nav is the number of pixel avalanches, cr G is the 
gain fluctuation parameter, <x el is the electronics noise inte- 
grated over At, and Run is the dark noise rate. Nav is re- 
lated to the number of photoelectrons at low light level by 
Nav = Npe(1 +Ncn), where Nqn the number of correlated noise 
avalanches per avalanche. This latter formula is an approxima- 
tion as it does not account for gain recovery and correlated noise 
avalanche created by other correlated noise avalanches. How- 
ever, the MC simulation includes both effects. Some conclu- 
sions can be drawn from this formula: 1) the integration gate 
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(Af) should be chosen so that Na v » Rdn^ in order to en- 
sure that dark noise does not contribute to the resolution, 2) the 
gain fluctuations do not contribute to the resolution significantly 
since ctq is only about 10%. This last conclusion highlights a 
significant difference between MPPCs and PMTs or standard 
Avalanche Photodiodes ( APDs), whose main contribution to the 
energy resolution arises from gain fluctuations. 

The simulated energy resolution is shown in Fig. 24 as 
a function of overvoltage with and without correlated noise 
(crosstalk and afterpulse). A gate of 540 ns was used to in- 
tegrate the charge. The light flash occurred 60 ns after the be- 
ginning of the gate and the photons were produced according 
to an exponential with a 7 ns time constant. The number of 
incident photons was set to 100 to match the average number 
of avalanches triggered by a minimum ionizing particle in T2K 
near detectors, which ranges between 20 and 35 avalanches. 
Without correlated noise the energy resolution would improve 
with increasing AV because of the increasing photodetection 
efficiency. In practice, when correlated noise is included the 
energy resolution reaches a minimum at AV = 1.8 V. Beyond 
1.8 V, correlated-noise-induced fluctuations worsen the energy 
resolution. Due to dynamic range constraints, in the T2K 
ND280 the MPPCs are operated at no more than AV = 1.33 V. 
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Figure 24: Simulated energy resolution as a function of overvoltage for a typical 
MIP signal of about 25 avalanches. The confounding effects include crosstalk, 
dark noise and afterpulses. The curve without these effects includes only the 
variation of the MPPC efficiency with overvoltage. 

The detector energy resolution is dominated by the photon 
counting statistics when AV is less than about 1.5 V; above 
1.5 V correlated noise contributes significantly. For photomulti- 
plier tubes and APDs, the contributions of gain fluctuations and 
correlated noise to the energy resolution are often assessed by 
calculating the excess noise factor (ENF). This better quantifies 
the contribution of the photosensor and the electronics system 
to the resolution by dividing out the fluctuations introduced by 
the photon counting statistics: 

F = oj f JN m (13) 

The dependence of the excess noise factor with overvoltage is 
shown in Fig. 25. The ENF increases with increasing overvolt- 
age following the increase of crosstalk and afterpulsing, which 
add additional avalanches in a stochastic manner. The ENF 



reaches 2 at a value of AV of about 1.5 V. The MPPC ENF 
is nevertheless significantly smaller than for APDs, whose ENF 
is always larger than 2 [20]. Overall, the MPPC contribution to 
the energy resolution is small for minimum ionizing particles 
that typically yield between 20 and 40 avalanches on average, 
even for T2K sub-detectors operating at AV = 1.33 V. 
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Figure 25: Excess Noise Factor as a function of overvoltage. 



4. Conclusion 

The T2K experiment ND280 complex of detectors uses a 
667-pixel MPPC developed by Hamamatsu Photonics specifi- 
cally for this experiment. It has a sensitive area of 1 .3x 1 .3 mm 2 
and a pixel size of 50x50 jj.m 2 ; the sensitive area is larger than 
those available previously and relaxes the mechanical toler- 
ances required for coupling to the WLS fibers used extensively 
in the experiment. We have performed detailed investigation 
of these devices and have developed an accurate model of the 
MPPC response to low light levels (where saturation effects can 
be neglected). 

MPPCs biased at the recommended Hamamatsu overvoltage 
(1.33 V) at T=25°C are characterized by the following parame- 
ters: photodetection efficiency of about 20% when illuminated 
with light from Yl 1 fibers (peaked in the green); a typical gain 
of 7.5x10 s ; the average dark rate is 700 kHz but can approach 
1 MHz; the crosstalk and afterpulse probability are estimated 
to be 9-12% and 14-16% respectively, with a combined total 
of 20-25%; and the recovery time constant of a single pixel 
is 13.4 ns. With such parameters, the device achieve the de- 
sired 20% energy resolution for for minimum ionizing parti- 
cles, which yield on average between 20 and 40 avalanches 
in the various components of the T2K near detector. Further- 
more, about 40,000 MPPCs were operated in the T2K neutrino 
beam in 2009-10 and no significant reliability issues were ex- 
perienced. 

Modeling the MPPC response by parameterizing dark noise, 
afterpulses, photodetection efficiency, crosstalk and gain varia- 
tion enables us to account for the contribution of the photosen- 
sor to the overall detector response accurately. The MPPC sat- 
uration effect should also be fully described by our simulations, 
but confirmation of this awaits additional controlled measure- 
ments for final validation. 
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